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ABSTRACT 

Several young supernova remnants exhibit thin X-ray bright rims of synchrotron radiation at their 
forward shocks. Thin rims require strong magnetic held amplihcation beyond simple shock compres¬ 
sion if rim widths are only limited by electron energy losses. But, magnetic held damping behind 
the shock could produce similarly thin rims with less extreme held amplihcation. Variation of rim 
width with energy may thus discriminate between competing inhuences on rim widths. We measured 
rim widths around Tycho’s supernova remnant in 5 energy bands using an archival 750 ks Chandra 
observation. Rims narrow with increasing energy and are well described by either loss-limited or 
damped scenarios, so X-ray rim width-energy dependence does not uniquely specify a model. But, 
radio counterparts to thin rims are not loss-limited and better rehect magnetic held structure. Joint 
radio and X-ray modeling favors magnetic damping in Tycho’s SNR with damping lengths ^1-5% of 
remnant radius and magnetic held strengths ^50-400 fiG assuming Bohm diffusion. X-ray rim widths 
are ~I% of remnant radius, somewhat smaller than inferred damping lengths. Electron energy losses 
are important in all models of X-ray rims, suggesting that the distinction between loss-limited and 
damped models is blurred in soft X-rays. All loss-limited and damping models require magnetic helds 
> 20 fiG, affirming the necessity of magnetic held amplihcation beyond simple compression. 

Subject headings: acceleration of particles — ISM: individual objects (Tycho’s SNR) — ISM: magnetic 
helds — ISM: supernova remnants — shock waves — X-rays: ISM 


1. INTRODUCTION 

Electrons accelerated in the forward shocks of young 
supernova remnants (SNRs) emit synchrotron radiation 
strongly in the shock’s immediate wake at radio wave¬ 
lengths and sometimes in X-rays. In a few cases, they 
quickly turn off downstream, producing a shell-like mor¬ 
phology of bright X-ra y and radio rims/hlaments due to 
line-o f-sight projection (Bamba et al. 2003 Reynoso et al. 
I997|) . Strong and timeTar i able synchrotron radiation 


e.g., Uchiyama et al. 2007 Patnaude & Eesen 2007), 


in conjunction with niu l tiwavelength spect r al modeling 


dAharonian et ah 2004 Acero et al. 2010 Ackermann 
et al. 2013), suggests that electrons are accelerated to 
i'eV energies in young SNRs. Although synchrotron 
emission due to accelerated electrons does not require ac¬ 
celeration of an unseen hadronic component, the prevail¬ 
ing theory of dihusive shock acceleration (DSA) should 
operate on both positive ions and electrons. Efficient 
hadron acceleration in supernova remnant shocks is a 
prime candidate source for galactic cosmic rays up to 
the co smic ray spectrum’s “knee” at around 3 PeV (Vink 
20121. But many fundamental questions about shock ac¬ 
celeration remain unanswered. Under what conditions 
do shocks accelerate particles efficie ntly? How are m ag- 
netic helds amplihed in such shocks? Reynolds (2008|) re¬ 
views relevant observations and open questions to date. 
These questions are relevant to ma ny astrophysical set - 


tings, such as Ear th’s bow shock (Elliso n et al. 1990) 

jets or active 


starburst galaxi es (Heckman et a l. 
galactic nuclei (Ichen et ai.||2014|), 


TMI) 


galaxy clusters (van 


_^_CRESST/Universit^_of^a2!t^2^j_College_Parkj_^D_207^ 

""^TRESST/Unlversitier^pa^ ResearchAss^iation 


Weeren et al.||2010|), and cosmological shocks (|Ryu et al. 

50DH|r -— - 

Spectral and spatial measurements of synchrotron rims 
can constrain downstream magnetic held strength and 
structure. If rim widths are set by electron energy losses, 
post-shock magnetic helds must be amplihed to qHOfjiG 
to account for the thinness of observed rims ( Vink fc 


Laming||2003[ IBamba et ^|2003H2005[JVblk et al. 

the 


Parizot et al.|2006[ ). In these models, the magnetic hel 

is assumed advected downstream and nearly constant 
over rim widths. Alternately, the magnetic held strength 
may be damped downstream of the shock and prevent 
electrons from radiating efficiently, so that thin rims re¬ 
hect magnetic held variation rather than e fficient par¬ 
ticle acceleration and synchrotron cooling (Pohl et al. 


20051. Damping, in particular, may permit less extreme 
magnetic held amplihcation. We refer to these as “loss- 
limited” and “damped” models for rim widths. We shall 
see that models can range continuously between these 
two cases, and that the distinction between the two can 
vary with observing frequency. 

The possi bility of damping in SNR sho cks has not been 
fully tested. Marcowith & Casse (2010) compared phys¬ 
ically motivated magnetic damping models to X-ray rim 
widths and synchrotron spectrum cut-offs and thus sug¬ 
gested that only young SNRs (age < 500 yr) can exhibit 
magnetic da mping in conjun c tion w ith efficient particle 


acceleration. Rettig & Pohl (2012) gave model predic¬ 


tions for several historical SNRs and proposed discrim¬ 
ination based on hlament spectra - the expectation is 
that damped spectra are softer, loss-limited harder. 

Hydrodynamic models can reproduce X-ray rim pro- 
hles reasonably well with both loss-limited and damped 
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magnetic field m odels ( QagsannQig^^ et al.|[2007 Mor- 


lino & Caprioli 2012 Slane et al. 2014). However, loss- 
limited models generally cannot reproduce thin radio 
rims. Radio-emitting GeV electrons do not lose substan¬ 
tial energy via radiation, so modeled intensities rise to a 
broad maximum toward the remnant interior, then grad¬ 
ually drop due to spheric i ty eff ects as the density drops 
in the interior. Reynolds (19881 empirically modeled ra¬ 
dio rims in the remnant oiSN 1006 and concluded that 
such thin rims required sharp gradients in electron energy 
density (from, e.g., time-v ariable particle accele r ation) 
or magnetic field strength. Cassam-Chenai' et al.| (20071 
used a 1-D hydrodynamic rnodel with nonlinear DSA to 
jointly model radio and X-ray rims in Tycho’s supernova 
remnant; neither loss-limited nor damped models could 
match radio profiles and radio/X-ray intensities simul¬ 
taneously. Moreover, gradual X-ray spectral variation 
observed downstream of th e shock front was poorly re - 
produced with both models ( Cassam-Chenai et al.|2007 ). 
But damped models were able generate limb- brightened 
radio rims at the forward shock, even if morphology was 
not en t irely c onsistent with observation. |Cassam-Chenai| 


et al. (2007) thus suggested that some combination of 
amplification and magnetic field variation might explain 
radio morp hology. 


Recently, Ressler et al. 


(|2014|) (hereafter, |R14| sought 
a rii 


to discriminate between damped and loss-lirnited rims by 
measuring rim width-energy dependence in X-ray ener¬ 
gies in SN 1006. In the simplest models, rim widths are 
expected to be roughly energy-independent if rims are 
damped, whereas widths should narrow with increasing 
energy if rims are energy loss-limited. R14 included a 
variety of effects which can blur t his di stinction. To fur¬ 
ther test these models, we follow |R14| by measuring X- 
ray rim widths at multiple energies in Tycho’s supernova 
remnant (hereafter, Tycho). Tycho exhibits an exten¬ 
sive shell of synchrotron-dominated thin rims around its 
periphery (Figure [^; the rims show very little thermal 
emiss ion, consistent with e xpansion into a low density 
ISM (|Williams et al.||20I3l). A deep 750 ks exposure of 
the entire remnant from 2009 allows fine sampling of the 
remnant rims. We build upon previous estimates of mag¬ 
netic field strength and particle diffusion in Tycho that 
draw from multiwavelength obse rvations and various as¬ 
sumptions on CR acceleration ( Volket ah 2002| 200^ 
Parizot et al.||200^ [Morlino &: Caprioli 2ol2[ ~ Rettig fc 


Pohl||2012[ e.g^ 


Uur procedure cl osely follows that of |R14[ We first 
review the model of |R14| used to model rim profiles and 
widths, then describe the procedure for selecting, mea¬ 
suring, and fitting rim widths to model width-energy de¬ 
pendence. We explore degeneracies in model fitting and 
consider radio rim morphology as an additional discrim¬ 
inant between models. 


2. NONTHERMAL RIM MODELING 
2.1. Particle transport 

The energy and space distribution of electrons at a 
supernova remnant’s forward shock controls the syn¬ 
chrotron rims we see in X-ray and radio. We assume 
that diffusive shock acceleration generates a power law 
distribution of electrons with an exponential cut-off at 
the forward shock and model 1-D steady-state plane ad- 



IC-Off 2ff 25”'0(r 0''24"'40‘ 

RA (J2000) 


Figure 1. RGB image of Tycho with region selections overlaid. 
Image bands are 0.7-1 keV (red), 1—2 keV (green) and 2-7 keV 
(blue). Bold region labels (1, 16) indicate region selections shown 
in Figures Filament 1: Regions 1—3, filament 2: regions 4-10, 
filament 3: regions 11-13, filament 4: regions 14-17, filament 5, 
regions 18—20. 


vection and diffusion of the electron distribution f{E^ x), 
where E is electron energy and x is distance downstream 
of the forward shock: 

in 


Following 

Berezhko & Volk 

2004); Cassam-Chenai et al. 

(2007 

|;1VJ 

iorlino et al. (2010 

TRettig & Fohll20l2|. The 


ing downstream of the shock; D is the diffusion coeffi¬ 
cient, Vd is fluid velocity downstream of the shock, and 
the constant b = 4e‘^/9my^ = 1.57 x 10“^ in appropri¬ 
ate CCS units arises from synchrotron power loss (i.e., 
dE/dt = —bB^E^, averaged over pitch angles). The ini¬ 
tial electron distribution is specified by an arbitrary nor- 
mali zatio n Kq, DSA cut-off energy Ecut (g iven in Sec¬ 
tion | 2.3 ), and spect ral index s = 2q; -|- 1. Zirakashvilij 


& Anaronian (2007) derive an electron energy spectrum 


with super-exponential cut-off , but we use a 

simp le exp onential cut-off for simplicity and consistency 
with R14| We have not yet specifie d th e cause of the 
injected spectrum cut-off (see Section 2.4), but the func¬ 


tional form of a power law with exponential cut-off is a 
good approximation to predictions for DSA spectra lim¬ 
ited by synchrotron losses, remnant age, or p article es¬ 


cape ( Webb et al.||1984 Reynolds 1998 2008). All con¬ 


stants and~equatIons are given in CCS (Gaussian) units. 
Our presentation is somewhat abbreviated, but a fuller 
exposition and literature review are given by |R14[ 

For Tycho , we a dopted radio spectral index a = 0.58 
(Sun et al. 2011) and hence electron spectral index 
s = 2a -I- 1 = 2.16. We assumed a remnant distance 


3 kpc (but cf. Hayato et al. 2010), which gives a shock 
radius of 1.08 x lO^ ’^ cm from the observed angular radius 
240" (|Green 2014) and further sets shock velocities for 
each rim pronle we consider. Tycho’s forward shock ve- 


locity varies with azimuth by up to a factor of 2 (Katsuda 
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et ah 2008); we linearly interpolated velocities reported 
^ williams et al. (2013) (rescaled to 3 kpc) to estimate 
individual shock velocities for each region. We assume a 
compression ratio of 4 as for a strong shock (unmodified 
by cosmic-ray pressure), and take downstream velocities 
Vd to be one-fourth the interpolated shock velocities. 

We assume isotropic diffusion and only consider par¬ 
ticle transport downstream of the forward shock. The 
velocity is assumed constant, as is the magnetic field 
for loss-limited model rims, in contrast to the expected 
Sedov-Taylor similarity solution for velocity and density 
in an adiabatic blastwave. Our assumptions of constant 
velocity, magnetic field, and plane flow should be reason¬ 
able as we generally consider synchrotron emission within 
10 % of the shock radius, r^, from the forward shock, 
though a few models are followed far enough inward that 
this approximation begins to break down. We consider 
profile emission strictly upstream o f both the contact dis - 
continuity and reverse shock (cf. Warren et al. 2005), 
and avoid regions where the contact discontinuity over¬ 
runs the forward shock. Our modeling neglects flow and 
electron spectrum modification due to the nearby con¬ 
tact discontinuity and particle acceleration at the reverse 
shock. These effects could matter (and we do not quan¬ 
tify their importance), but we expect that at Tycho’s age. 
X-ray emission is dominated by forward shock transport. 
Given the uncertainty in shock morphology, our model 
seeks only to capture the most relevant physics. More 
sophisticated work may treat, e.g., sphericity, shock pre¬ 
cursors, anisotr opic diffusion, and in i ection/acceleration 


efficiency (e.g. , Rev ille & Bell 2013 Bykov et al. 2014 
Ferrand et al.|[20T4 and references therein). 


To determihe rim profiles and widths, we compute the 
electron distribution using G reen’ s function solut i ons by 
Lerche & Schlickeiser (1980) and Rettig & Pohl (2012), 
with the caveat that D{x)B’^{x) is assumed constant; 
we discuss this assu mption further below. The solutions 
are fully described in |R14| using notation similar to ours. 
The electron distribution may be integrated over the one- 
particle synchrotron emissivity G{y) to obtain the “to¬ 
tal” emissivity: 


jv{x)(x / G{y)f{E,x)dE 


( 2 ) 


where y = u/{ciE^B) is a scaled synchrotron frequency 
and G{y) = y K^/^{z)dz with K^/^{z) a modified 
Bessel function of the second kind ( jPacholczyk 1970); 
the constant ci = 6.27 x 10^® in GGS units. Integrat¬ 
ing emissivity over lines of sight for a spherical remnant 
yields intensity as a function of radial coordinate r: 


I„{r) = 2 j ‘ > (ts - \/s2~+r2^ 


ds 


( 3 ) 


where s is the line-of-sight coordinate and is shock ra¬ 
dius. We take the full width at half maximum (FWHM) 
of the resulting intensity profile as our metric for modeled 
rim widths. Using FWHM as opposed to, e.g., full width 
at three-quarters maximum, excludes measured rims and 
model parameters where X-ray intensity does not drop to 
half maximum immediately behind the rim. Our results 
thus focus on the most well-defined rims in Tycho rather 
than the global shock structure. 


2.2. Magnetic fields and damping 

We consider two scenarios for post-shock magnetic 
field: (1) a constant field B{x) = Bq corresponding to 
loss-limited rims, and ( 2 ) an exponentially damped field 
of form: 

B{x) = {Bo - Hmin) exp {-x/ab) + Hmin (4) 

following ( Pohl et al.|2005 ). Here Bq is the magnetic field 
immediately downstream of the shock, i.e. Bq = B{x = 
0 ), and Ub is an e-folding damping lengthscale; our use of 
Bq for downstream magnetic field departs from typical 
notation. A typic al lengthscale for ab is 10^® to 10^^ cm 


( Pohl et ^ 12665 ), corresponding to ~0.1-1% of Tycho’s 
radius. Hereafter, we report ab in units of shock radius 
Tg unless otherwise stated. 

The distinction between damped and loss-limited rims 
is somewhat arbitrary; as Of, — 00 , model results con¬ 
verge to loss-limited rims. Moreover, rims much thin¬ 
ner than the damping length are effectively loss-limited 
as electrons radiate in a nearly constant magnetic field. 
Furthermore, the dependence on observing frequency of 
electron losses and diffusion means that a model may be 
damping limited at one frequency and loss-limited at an¬ 
other. In the following analysis, we deem all fits with 
finite ab to be damped, but we compare rim widths and 
damping lengths from such fits further below to better 
distinguish damped and loss-limited rim behavior, set by 
a combination of ab, Bq, and other model parameters. 

2.3. Diffusion coefficient 

Most previous work has assumed Bohm-like diffusion 
in plasma downstream of SNR shocks. Bohm diffusion 
assumes that the particle mean free path A is equal to 
the gyroradius Vg = E/{eB), yielding diffusion coefficient 
= Ac/3 = cE/{3eB)-, here c is the speed of light, E 
is particle energy, e is the elementary charge, and B is 
magnetic field. Bohm-like diffusion encapsulates diffu¬ 
sion scalings oi D (x E, introducing a free prefactor 77 
such that A = rjCg allows for varying diffusion strength. 
However, Bohm diffusion at 77 = 1 is commonly con¬ 
sidered a lower limit on the diffusion coefficient at all 
energies. 

We consider a generalized diffusion coefficient with ar- 
bitr ary power law depe ndence upon energy following, 
e.g., Parizot et al. (2006): 


D{E) = ^ = 77,Ub {Eb) 


E 


( 5 ) 


where /i parameterizes diffusion-energy scaling and 77 now 
has units of erg^”'^. The right-hand side of equation (1^ 
introduces rjh, a dimensionless diffusion coefficient scalra 
to the Bohm value at a fiducial particle energy Eh- Note 
that rjh and 77 are related as 77 = r]fi{Eh)^~^, and rj = rjh 
for Bohm-like diffusion (/r = 1). For subsequent analy¬ 
sis, we take fiducial electron energy E 2 = E^ correspond¬ 
ing to a 2 keV synchrotron photon and report results in 
terms of 772 = rjh- Although E 2 varies with magnetic field 
as E 2 oc and thus 772 may vary around Tycho’s 

shock for /r 7 ^ 1 , tying rjh to a fixed observation energy 
gives a convenient sense of diffusion strength regardless 
of the underlying electron energies. 

The solutions t o equation (f^ given by Lerche & 
Schlickeiser (1980) assume D{x)B'^{x) constant to ren- 
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der equation Q semi-analytically tractable. Although 
this assumption has no obvious physical basis, it con¬ 
tains the qualitatively correct behavior D constant if B 
is constant and D smaller for larger B. We enforce it by 
modifying th e diffusion coefScient in the damping model 
as, following Rettig & Pohl (2012|: 


D{E,x) = 


rjcE^ 

SeBn 


Brai 


Bq — i?i 
Bn 


min^-x/ab 


-2 


( 6 ) 


This strengthens the spatial-dependence of the diffusion 
coefficient as compared to the expected D{x) oc l/i?(a;). 

2.4. Electron energy cut-off 

We assume that the DSA process is limited by syn¬ 
chrotron losses at high energies and hence determine 
the Ecut by equating synchrotron loss and diffusive ac¬ 
celeration timescales. Here we ignore the pos sibility of 
age-li mited or escape-limited acceleration (e.g.,|Reynolds' 
19981. In the former case, low magnetic-field strengths 
could mean that Tycho’s age is less than a synchrotron 
loss time; the maximum energy is obtained by equating 
the remnant age and acceleration timescale. In the lat¬ 
ter case, the diffusion coefficient upstream may increase 
substantially above some electron energy due to an ab¬ 
sence of appropriate MHD waves. However, the mag¬ 
netic field strengths we find below justify the assumption 
of loss-limited acceleration. For low energies and small 
synchrotron losses (cooling time longer than acceleration 
time), electrons are efficiently accelerated; near or above 
the cut-off energy, electrons will radiate or escape too 
rapidly to be accelerated to higher energies and the en¬ 
ergy spectrum drops off steeply. The cut-off energy is 
given as: 


observable to estimate shock diffusion and is given by: 
J^cut = Cjn (13.3 erg)!^ (100 /rG) (2657 erg^) 


10 ® cm/s 


im 


1 + fi 


( 8 ) 


The cut-off frequency is independent of magnetic field 
B for all values of /i (but recall that the electron en- 
ergies associat e d wit h 772 will depend on B for 77 7 ^ 1 ). 
Parizot et al. (2006| previously used measurements of 
synchrotron cut-otts to estimate diffusion coefficients in 
Tycho and other historical supernova remnants. 

We point out that the electron spectrum softens down¬ 
stream of the shock due to synchrotron losses, but the 
local spectrum at any radial position will be a ste eply 
cut-off power law ( Webb et al.|1984 Reynolds|1998 |. No 
steepening from to is observed because the 

cut-off limits the electron spectrum at high energies. A 
homogeneous source of age t in which electrons are con¬ 
tinuously accelerated throughout, with an initial straight 
power law distribution to infinite energy, will produce a 
steepened power law distribution above the energy at 
which the synchrotron loss time equals the acceleration 
time (Kardashev 19621. If we model emission without 
an initial exponential cut-off (i.e., inject a straight power 
law) in a constant magnetic field, then integrated spectra 
in our model would steepen by about one power. But, 
those assumptions do not apply to the current situation 
of continuous advection of electrons, with an energy dis¬ 
tribution which is already a cut-off power law, through 
a regi on of non-constant magnetic field. See [Reynolds] 
(20091 for a fuller discussion of synchrotron losses in non- 


homogeneous sources. 


Beut = (8.3 TeV)"/(^+^) 


Bn 


( 


100 fiG 
\2/(i-eM) 




10 ® cm s“i / 


V 




(7) 


This result is derived by Parizot et al. (2006) for /r = 1 
assuming a strong shock with compression ratio 4 and 
isotropic magnetic turbulence both upstream and down¬ 
stream of the shock. 

The cut-off energy and accelerated electron spectrum 
depend on the magnetic field, which varies in a damped 
model. Marcowith & Casse (2010) show that various 
turbulent damping mechanisms can modify the acceler¬ 
ated spectrum, as particles traveling downstream may 
not be effectively reflected back across the shock and 
further accelerated. Nevertheless, we make the simpli¬ 
fying assumption that particle acceleration is controlled 
by diffusion at the shock and neglect spatially varying 
diffusion and magnetic fields in the acceleration process; 
equation 0 stands as evaluated with shock magnetic 
field strength Bq. Cut-off energies of ~1-10 TeV can 
be plausibly achieved in the presenc e of damping due to 
Alfve n and magneto-sonic cascades (Marcowith & Casse 

Mo|. 


T^s the DSA imposed electron cut-off results in a cut¬ 
off of SNR synchrotron flux, the synchrotron cut-off fre- 
que ncy i'„,t = CmE^^t B with Cm = 1.82 x 10 ^® in cgs units 
(e.g. Pacholcz'^ 1970), which is the peak frequency emit¬ 
ted by electrons of energy Bcut, provides an independent 


2.5. Rim width-energy dependence 

X-ray rim widths are controlled by synchrotron losses, 
particle transport, and magnetic fields immediately 
downstream of the shock, each of which influences rim 
width-energy scaling differently. If diffusion is negligi¬ 
ble and the downstream magnetic field is constant, loss- 
limited rims narrow with increasing energy as more ener¬ 
getic electrons radiate and cool more quickly. But diffu¬ 
sion will dilute this effect: more energetic electrons may 
diffuse further upstream or downstream than would be 
expected from pure advection, smearing out rims and 
weakening energy dependence at higher energies. Mag¬ 
netic fields damped on a length scale comparable to hla- 
ment widths should also weaken rim width-energy depen¬ 
dence - if the magnetic field turns off, synchrotron radia¬ 
tion turns off regardless of electron energy. Additionally, 
once Bq varies, one observes electrons of different energy 
at different distances behind the shock. 

Figure plots model X-ray and radio profiles for a 
range of Bq and at, values to illustrate how damping and 
magnetic field strength impact width-energy dependence, 
which we will now explore. We discuss and incorporate 
model radio profiles into our analysis in Section 

2.5.1. Undamped models 

Following |R14[ we parameterize rim width-energy de¬ 
pendence in terms of a scaling exponent toe defined as: 

w{iy) oc t/"*® (9) 
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So = 20 [iG 


Bo = 30 |lG 


Bo = 100|iG 



Radial position (arcsec.) 



Radial position (aicsec.) 


Radial position (arcsec.) 


Radial position (arcsec.) 


Figure 2. Radio and X-ray rim profiles for a range of magnetic fields (Bq) and damping lengths (a^, shaded regions) with /i = 1, 772 = 1, 
and -Bjnin = 5 mG fixed. X-ray rim energies (0.7, 2, 4 keV) are representative of energy bands used in our rim width measurements. This 
plot summarizes several key features of our model: (1) radio profiles (dashed red) are strongly affected by damping for all values of Bq 
because synchrotron losses are negligible at low electron energies. (2) X-ray profiles (solid green, black, blue) are influenced by synchrotron 
losses even in the presence of strong damping, which can be seen going from left to right (increasing Bq). (3) When X-ray rim widths are 
smaller than field damping only weakly affects rim widths and width-energy dependence (top right panels). (4) Strongly damped X-ray 
rims can show significant width-energy dependence in our model due to synchrotron losses beyond a^,. This occurs for smaller Bq where 
the contrast between Bq and Smin is less extreme (bottom left panels). 


where w{v) is filament FWHM as a function of observed 
photon frequency v, and the exponent ttie = toe(^') is 
energy dependent. We may refer to observed photons 
interchangeably by energy or frequency v, but E is re¬ 
served for electron energy. 

To better intuit the effects of advection and diffu¬ 
sion on rim widths, we introduce advective and diffusive 
lengthscales for bulk electron transport. These depend 
on electron energy; we write them in terms of the peak 
frequency radiated by electrons of energy E,v = CmE^B. 

kd = VdTsynch oc (10) 

^diff = VT'7'synch oc 7 yl/ 25 -(/"+ 5 )/ 4 ^(M-l )/4 

The characteristic time is the synchrotron cooling time 
Tsynch = ll{bB^E) with b = 1.57 X 10“^. For /r = 1, /diff 

is independent of v and both Zdiff and kd scale as Bq . 
If both diffusion and magnetic field damping are negli¬ 
gible and electrons are only loss-limited as they advect 


downstream, attains a minimum value ttie = —1/2 
as rim widths are set by kd > At higher energies 
where Zdiff > ^adj diffusion increases toe from —1/2 to a 
value between —1/4 and 1/4 for ^ = 0 and 2 respectively 
(R14 Figure 3). The presence of an electron energy cut- 
ott decreases toe slightly in all cases due to the decreased 
numbe r of e lectrons and hence thinner rims at higher en¬ 
ergies (R14 Figure 5), but the qualitative behavior is the 


same. 


2.5.2. Field damping effects 

We expect magnetic damping to produce compara¬ 
tively energy-independent rim widths. If synchrotron 
rim widths are set by magnetic damping at some ob¬ 
servation energy, then rims will be damped at all lower 
observation energies as well. Then rim widths will be 
relatively constant (small |toe|) below a threshold en¬ 
ergy and may decrease, or even increase once advection 
and/or diffusion control rim widths at higher photon en- 
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ergies (advection: 4d < cLb] diffusion: Zdiff > la.d,a-b)- 
Thus, we intuit that rim widths should roughly scale 
as w ~ min (of), max (/ad, ^diff))- This is correct except 
for one key region of parameter space: strong damp¬ 
ing with weak magnetic held amplihcation, where syn¬ 
chrotron losses downstream of the FWHM create energy 
dependent widths even when at ^ /ad, ^diff- 
This counter-intuitive energy dependence for strongly 
damped models occurs when rim brightness remains 
above half-maximum within ~a{, of the shock. Farther 
downstream, synchrotron losses in the reduced magnetic 
held drive intensity down to half-maximum; losses 

in a nearly constant held cause energy dependence de¬ 
spite magnetic damping at the shock. Synchrotron losses 
far from the shock will only affect the spectrum at higher 
frequencies; at lower frequencies, the spectrum will be 
effectively constant after the magnetic field’s initial de¬ 
cay. Thus, there will be a characteristic frequency be¬ 
low which the FWFIM ceases to be dehned. The energy 
dependence of rim FWFIMs becomes large and diverges 
near this frequency. The characteristic frequency and 
value of Toe are sensitive to our dehnition of rim width 
(FWHM), but the physical behavior we describe (diver¬ 
gence of ruE at characteristic frequency, for appropriate 
damping parameters) will occur regardless of our dehni¬ 
tion of rim width. 


3. OBSERVATIONS 
3.1. Data and region selections 

We measured synchrotron rim full widths at half max¬ 
imum (FWHMs) from an archival Chandra ACIS-I ob¬ 
servation of Tycho (RA: 00*'25“19!0, dec: -k64°08'10"0; 
J2000) between 2009 Apr 11 and 2009 May 5 (PI: J. 
Hugh e s: Ob slDs: 10093-10097, 10902-10906); [Eriksen 
et al. (20111 present additional observation information. 
The total exposure time was 734 ks. Level 1 Chan¬ 
dra data were reprocessed with CIAO 4.6 and CALDB 
4.6.1.1 and kept unbinned with ACIS spatial resolution 
0.492". Merged and corrected events were divided into 
hve energy bands: 0.7-1 keV, 1-1.7 keV, 2-3 keV, 3-4.5 
keV, and 4.5-7 keV. We excluded the 1.7-2 keV energy 
range to avoid Si XIII (Hea) emission prevalent in the 
remnant’s thermal ejecta which might contaminate our 
nonthermal prohle measurements. 

We selected 20 regions for prohle extraction around 
Tycho’s shock (Figure based on the following criteria: 
(1) hlaments should be clear of spatial plumes of thermal 
ejecta in Chandra images, which rules out, e.g., areas of 
strong thermal emission on Tycho’s eastern limb; (2) hl¬ 
aments should be singular and localized, so multiple hla¬ 
ments should either not overlap or completely overlap; 
(3) hlament peaks should be evident above the back¬ 
ground signal or downstream thermal emission (rules 
out faint southern hlaments). We accepted several re¬ 
gions with poor quality peaks in the lowest energy band 
(0.7-1 keV) so long as peaks in all higher energy bands 
were clear and well-ht. We grouped regions into 5 hla¬ 
ments by visual inspection of the remnant. Within each 
hlament, we chose region widths to obtain comparable 
counts at the thin rim peak. All measured rim widths are 
at least 1". The narrowest rim widths may be slightly 
over-estimated due to the Chandra point-spread func¬ 
tion (PSF) at ^4' off the optical axis, which has FWHM 


< 1.4" at 6.4 keV and < 1" at 1.5 keV. For simplicity, 
we neglect PSF ehects in our analysis. Our observations 
of rim width-energy dependence are thus somewhat con¬ 
servative at the highest energy 


3.2. Filament spectra 


We extracted spectra at and immediately behind thin 
rims in each region (“rim”, “downstream” spectra respec¬ 
tively) to conhrm that rim width measurements are not 
contaminated by thermal line emission. The two extrac¬ 
tion regions are deter min ed by our empirical fits of rim 
profile shape (Section |3.3| ). The rim section is the small¬ 
est sub-region containing the measured FWHM bounds 
from all energy bands. The downstream section extends 
from the interior thin rim FWHM to the intensity mini¬ 
mum behind the rim (specifically, the do wnst ream profile 
fit domain bound described in Section 3.3). To illus¬ 
trate our selections. Figure plots example rim profiles 
(4.5-7 keV) with the downstream and rim sections high¬ 
lighted. 

Spectra were binned to a minimum of 15 cts/bin. We 
extracted background spectra from circular regions (ra¬ 
dius ~30") around the remnant’s exterior; for each re¬ 
gion’s rim and downstream spectra, we subtracted the 
closest background region’s spectrum. 

We fit each region’s rim and downstream spectra 
to an absorbed power law model (XSPEC 12.8.1, 
phabs*powerlaw) between 0.5-7 keV with photon index 
F, hydrogen column density Nn, and a normalization as 
free parameters. Tablelists best fit parameters and re¬ 
duced values for all regions. Rim spectra are well-fit 
by the power law model alone; the best fit photon in¬ 
dices (2.4-3) and column densities (0.6-0.8 x 10^^ cm“^) 
are consistent with previous spectral fits to Tycho’s non- 


thermal rims (Hwang et al. 2002[ Cassam-Chenai et al. 


20071. 

Downstream spectra are poorly ht by the absorbed 
power law model due to thermal contamination from 
Si XIII and S XV Hea line emission at 1.85 and 2.45 
keV. To conhrm that thermal emission is dominated by 
these two lines near the shock, we also performed hts with 
(1) both lines excised (1.7-2.0 keV, 2.3-2.6 keV counts 
removed) and (2) with both lines htted to Gaussian pro- 


hles. Fits with lines excised yield Xred values between 
1-5. Fits with lines htted to Gaussian prohles yield x^ed 
values 0.83-1.6. In both hts (lines excised or modeled), 
we hnd somewhat smaller best ht column densities (0.3- 
0.8 X 10^^ cm“^) but similar best ht photon indices (2.6- 
3.1), compared to those of the rim spectra. The consis¬ 
tent photon indices indicate that the same synchrotron 
continuum is present beneath thermal li ne em ission. 

We also htted “rim” spectra (Section |3.2[ ) to the ab¬ 
sorbed XSPEG model srcut, modihed to ht in log- 
frequency space, srcut models a power law X-ray syn¬ 
chrotron spectrum set by a radio spectral index a with 



lists best spectrum ht parameters for each region. The 
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Figure 3. Spectra and fits from Regions 1 (top) and 16 (bottom) show varying rim morphology; Region 1 shows a rim where the 0.7-1 keV 
peak could not be fit. Left: 4.5-7 keV profiles with downstream (blue) and rim (grey) sections highlighted. Intensity is in arbitrary units 
(a.u.). Middle: downstream spectra with absorbed power law fit; Si and S lines at 1.85, 2.45 keV are clearly visible. Right: rim spectra 
with absorbed power law fit show that rims in each region are likely free of thermal line emission. 


Table 1 

Absorbed power law spectrum fit parameters 


Downstream, power-law Rim, power-law Rim, srcut 


Region 

Nh 

r 

X?ed (dof) 

Nh 

r 

X?ed (dof) 

Nh 

I^CUt 

X?ed (dof) 



(-) 


(lO^^cm”^) 

(-) 


(10^^cm“^) 

(keV/fi) 


1 

0.53 

2.72 

2.25 (186) 

0.72 

2.84 

1.20 (284) 

0.62 

0.30 

1.17 (284) 

2 

0.68 

2.99 

5.07 (178) 

0.69 

2.83 

1.10 (202) 

0.60 

0.30 

1.08 (202) 

3 

0.67 

2.96 

1.97 (186) 

0.77 

2.80 

1.15 (167) 

0.68 

0.33 

1.14 (167) 

4 

0.59 

2.97 

1.43 (163) 

0.70 

2.84 

1.21 (278) 

0.61 

0.29 

1.15 (278) 

5 

0.62 

2.93 

4.76 (265) 

0.73 

2.88 

1.18 (255) 

0.64 

0.27 

1.11 (255) 

6 

0.68 

3.00 

2.12 (200) 

0.74 

2.85 

0.96 (231) 

0.64 

0.29 

0.93 (231) 

7 

0.65 

3.02 

1.00 (142) 

0.82 

2.97 

1.14 (224) 

0.71 

0.23 

1.14 (224) 

8 

0.74 

2.93 

1.38 (170) 

0.75 

2.71 

0.98 (198) 

0.66 

0.41 

0.96 (198) 

9 

0.78 

3.03 

1.12 (157) 

0.82 

2.83 

0.90 (175) 

0.73 

0.30 

0.88 (175) 

10 

0.62 

2.86 

1.40 (220) 

0.77 

2.76 

0.98 (164) 

0.68 

0.36 

0.96 (164) 

11 

0.67 

2.94 

2.56 (137) 

0.69 

2.60 

1.10 (153) 

0.61 

0.55 

1.07 (153) 

12 

0.61 

2.79 

2.65 (137) 

0.64 

2.44 

0.90 (172) 

0.57 

0.88 

0.90 (172) 

13 

0.61 

2.98 

3.12 (198) 

0.67 

2.73 

1.12 (235) 

0.59 

0.38 

1.09 (235) 

14 

0.46 

2.93 

1.37 (148) 

0.63 

2.93 

0.96 (167) 

0.54 

0.23 

0.95 (167) 

15 

0.43 

2.92 

1.33 (150) 

0.65 

2.84 

1.05 (183) 

0.57 

0.29 

1.03 (183) 

16 

0.48 

2.94 

2.04 (189) 

0.67 

2.80 

1.13 (182) 

0.58 

0.32 

1.12 (182) 

17 

0.48 

2.86 

1.70 (188) 

0.68 

2.83 

0.97 (187) 

0.59 

0.30 

0.94 (187) 

18 

0.44 

2.87 

1.91 (200) 

0.64 

3.02 

1.20 (220) 

0.55 

0.19 

1.13 (220) 

19 

0.40 

2.84 

1.31 (133) 

0.66 

2.78 

1.01 (157) 

0.57 

0.34 

0.99 (157) 

20 

0.40 

2.75 

3.01 (140) 

0.63 

2.81 

1.11 (192) 

0.55 

0.31 

1.07 (192) 


Note. — Absorbed power law fit parameters are photon index F and hydrogen column density Nu- srcut fits performed in log-frequency space; 
h is Planck’s constant and Vcut ^ cut-off frequency. Horizontal rules group individual regions into filaments. 



















































Tran et al. (September 4, 2015) 


fitted cut-o ff frequency is typic ally 0.3 keV/h, consistent 
with fits by Hwang et al. (2002), but Regions 11, 12 have 
unusually high cut-oft frequencies 0.55, 0.88 keV/h con¬ 
sistent with harder rim spectra. 

Our spectral fitting confirms that all selected region 
are practically free of thermal line emission, as already 
suggested by visual inspection (Figurej^. Excluding 1.7- 
2 keV photons in rim width measurements further lim¬ 
its thermal contamination as 1.85 keV Si line emission 
is over a third of Tycho’s th ermal flux as detected by 
Chandra ( Hwang et al.|[2002 |. 


3.3. Filament width measurements 

We obtained radial intensity profiles in five energy 
bands from ~10-20" behind the shock to ~5-10" in front 
for each region. To increase signal-to-noise, we integrate 
along the shock (5-23") in each region. Plotted and 
fitted profiles are reported in vignetting and exposure- 
corrected intensity units; error bars were computed from 
raw counts assuming Poisson statistics. Intensity profiles 
peak sharply within ~2-3" behind the shock, demarcat¬ 
ing the thin rims, then fall off gradually until thermal 
emission picks up further behind the shock. 

We fitted rim profiles to a piecewise two-exponential 
model: 

/i(r) = (12) 

|Hdexp(^J+Q, r<ro 

where h{r) is profile height and r is radial distance from 
remnant center. The rim model has 6 free parameters: 

^d-> ^u-> iind C(i^ A(i {Cii Cd^ en- 

fo rces continuity at r = rn. Our model is similar to that 
of lBarn ba et al. (2003 20051 and differs slightly from that 
of |R14 To ht only the nonthermal rim in each intensity 
profile, we selected the fit domain for each profile as fol¬ 
lows. The downstream bound was set at the first local 
data minimum downstream of the rim peak, identified by 
smoothing the profiles with a 21-point (^10") Hanning 
window. The upstream bound was set at the profile’s 
outer edge. Figure [^illustrates the fit domain selections 
for two example regions. 

From the fitted profiles we extracted a full width at 
half maximum (FWHM) for each region and each en¬ 
ergy band after subtracting a constant background term 
min(C'u, Cd). We could not resolve a FWHM in 8 of 20 
regions at 0.7-1 keV (Table [2); in these regions, either 
the downstream FWHM bourn would extend outside the 
fit d oma in or we could not find an acceptable fit to equa¬ 
tion ( |l^ . We were able to resolve FWHMs for all regions 
at higlier energy bands (1-7 keV). 

To estimate FWHM uncertainties, we horizontally 
stretched each best fit profile by mapping radial coor¬ 
dinate r to r'(r) = r(l -|- ^(r — ro)/(50" — tq)) with 
^ an arbitrary stretch par ame ter and rg the best fit 
rim center from equation (12), yielding a new profile 
h'{r) = h{r'{r)). We varied ^T^nd hence rim FWHM) to 
vary each profile fit by 2.7 and took stretched FWHMs 
as upper/lower bounds on reported FWHMs. This pro¬ 
cedure again follows |RT^ 

3.4. Filament model fitting 


We fit model FWHM predictions given by equation ([^ 
to measured rim widths as a function of energy by vary¬ 
ing several physical parameters: magnetic field strength 
Bq, normalized diffusion coefficient rj 2 , diffusion-energy 
scaling exponent /x, and minimum field strength Hmin 
and lengthscale ab for a damped magnetic field. We 
mapped each width measurement to the lower energy 
limit of its energy band; e.g., 0.7-1 keV is assigned to 
0.7 keV and fitted to model profile widths at 0.7 keV. 
Width errors in our least squares fits average the positive 
and negative errors on each FWHM measurement. For a 
given set of model parameters, we numerically computed 
intensit y pr ofiles and hence model FWHMs as detailed in 
Section [TT| we then used a Levenberg-Marquardt fitter 
to seek model parameters yielding best fit FWHMs. To 
assist the nonlinear fitting, we tabulated model FWHM 
values on a large grid of model parameters and used best 
fit grid parameters as initial guesses for fitting. We re¬ 
quired 772 to be positive and deemed best fit values with 
772 > 10® and Bq > 10 mG to be effectively uncon¬ 
strained. In subsequent analysis we focus on fits with 
/X = I and 772 = 1 fixed, though we discuss the effects of 
varying both 77 and 772 (and fitting for 772 ) as well (recall 
that for /X = 1, ?7 is energy-independent). 

The purely loss-limited model (constant downstream 
magnetic field Bq) has three parameters /x, 772 , and Bq. 
To make nonlinear fitting tractable, we fixed /x in all 
fits and considered /x = 0, 1/3, 1/2, 1, 1.5, and 2. In 
particular, nonlinear diffusion-energy scalings with 77 = 
1/3 and 1/2 may arise from Kolmogor ov and Kraichn an 
turbulent energy spectra respectively ( Reynolds||2004 ). 

For damped magnetic field fits, we held the remnant 
interior field strength Rmin constant at 5 txG, sl ightly 
higher than typical intergalactic values of ~2-3 txG (jLyne 
'& Smith||1989 |Han et al.||2006|). We stepped at through 
14 different values between U.5 and 0.002 (sampling most 
finely between 0.01 and 0.002). To ensure that damped 
fits generate rims influenced by magnetic damping, we ar¬ 
bitrarily require that the rim FWHM at 2 keV be strictly 
greater than the fit value of at,; we revisit this require¬ 
ment below. For best fits, we report the value of at, 
yielding the smallest value, with the caveat that that 
our at, sampling is relatively coarse. 

Predicted rim widths are subject to resolution error 
in the numerical integrals (discretization over radial co¬ 
ordinate, line-of-sight coordinate, electron distribution. 
Green’s function integrals). We chose integration reso¬ 
lutions such that the fractional error in model FWHMs 
associated with halving or doubling each integration res¬ 
olution is less than 1% for the parameter space relevant to 
our filaments. The maximum resolution errors in a sam¬ 
ple of parameter space are typically 0.1-1%, but mean 
and median errors are typically an order of magnitude 
smaller than maximum errors. 


4. RESULTS 
4.1. Rim widths 

Measured rim widths decrease with energy in most re¬ 
gions and energy bands. Table reports FWHM mea¬ 
surements for all of our regions. We also report ttxe values 
for all but the lowest energy band, computed point-to- 
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Figure 4. Best fit profiles with measured FWHMs demarcated for each energy band in Region 1 (top) and Region 16 (bottom). Hereafter, 
we use Regions 1 and 16 to illustrate results for profiles of differing quality and absolute width. We could not measure a 0.7-1 keV FWHM 
in Region 1, reflected in TableData in red were excluded from profile fit domains as described in text. 


point between discrete energy bands as: 


mE{E2) 


\w{w2/wi) 

HE2/EI) 


(13) 


where wi,W 2 and Ei,E 2 are FWHMs and lower energy 
values for each energy band - e.g., ttie at 1 keV is com¬ 
puted using FWHMs from 0.7-1 keV and 1-1.7 keV, 
with i ?2 = 1 keV and E^ = 0.7 keV. Errors on ttie are 
propagated in quadrature from adjacent FWHM mea¬ 
surements. 

Although the measurement scatter is quite large, 
shown dramatically in the point-wise computed toe val¬ 
ues, the mean rim width decreases consistently with in¬ 
creasing energy. Furthermore, mean rriE values are con¬ 
sistently negative and tend smoothly towards 0 (weaker 
energy-dependence) with increasing energy. Errors on 
FWHM measurements are typically < 10%, reflecting 
the high quality of the underlying Chandra data. Scat¬ 
ter in FWHM measurements may be attributed in part 
to ( 1 ) our measurement procedure, which depends on an 
empirical choice of profile fit function, and ( 2 ) variation 
in Tycho’s rim morphology (e.g.. Figure^. 

4.2. Model fit results 

Best fit results for loss-limited and damped cases are 
summarized in Figure and Table with fj, = 1 and 
772 = 1 both fixed. In ^ tables, we report values of 
rather than Xred comparison of different fits, as 

we often manually stepped parameters that were held 
constant in fitting. 

Table reports model fits with 772 fixed to values de¬ 
rived from srcut spectrum fits for i/cut and hence 772 
from equation ®, with Bq varying freely; cutoff-derived 
772 values of ^ITT require increased magnetic fields com¬ 
pared to the 772 = 1 case if rims are purely loss-limited. 
We derive different 772 values for varying 77 and find that 
fitted Bq varies only weakly with fj,] for example, Re¬ 
gions 1 and 16 have best fit values of Bq increasing over 
252-266 /rG and 660-700 fiG respectively as 77 ranges 
between 0 - 2 , for loss-limited fits. 

Magnetically damped rims are able to fit width-energy 
dependence in our measurements at least as well as 
purely loss-limited rims for an acceptable range of at 
values. As loss-limited models are a subset of damping 


models, this is expected. But fits with damping lengths 
small enough to exert influence on rim widths are also 
permissible. In both loss-limited and damped fits, we do 
not observe systematic variation of fit parameters with 
azimuth around the remnant. 

4.3. Fitting parameter degeneracy 

We briefly explore the effects of varying 772 and other fit 
parameters, before focusing on results with 77 = 772 = 1 
for simplicity. Fixing 772 from measured i/cut values does 
not yield any obvious insight. But, diffusion coefficients 
computed in this manner may give more credible derived 
estimates of Bq if the DSA assumptions invoked are ac¬ 
curate. 

Loss-limited model values of Bq and 772 may covary 
without strongly altering fit quality in many regions. If 
diffusion is the primary control on rim width - i.e., /diff > 
^ad “ then 772 and Bq become degenerate as the product 


-3/2 


(equation (111). If advection is the primary control on 


exerts most control on rim width w ~ Zdiff 


rim width (widths narrow rapidly with energy; i.e., rriE 
—0.5), then 772 I becomes unimportant and fits are 
well-behaved with effectively one free parameter. 

We also performed loss-limited model fits with 77 fixed 
between 0 and 2 and both 772 and Bq free. Fits with 
77 < 1 generally yield larger parameter values and errors 
for both 772 and Bq, and they are more likely to be ill- 
constrained (e.g., 772 > 10 ^ or Bq > 10 ^ fiG). Fits to the 
same data with varying 77 can yield 772 varying by 1-2 or¬ 
ders of magnitude. The x^ values for individual fits are 
variable and large (^ 1 ), so we cannot favor or disfavor 
particular values of 77 and 772 . But, neglecting the mag¬ 
nitude of our x^ values, values of 77 > 1 are qualitatively 
favored by x^ in niost regions. This trend may be par¬ 
tially an artifact of the correlation between Bq and 772 as 
they are not entirely independent parameters, but they 
are less correlated for larger 77 ; fits at smaller 77 may have 
fewer (non-integer) degrees of freedom, partially offset¬ 
ting our observations. 

Damping modifies the degeneracy in Bq and 772 . For 
small values of ab (strong damping), increased diffu¬ 
sion ( 772 ) can cause rim widths at all energies to narrow 
counterintuitively. Speculatively, spatial variation of the 
diffusion coefficient may oppose downstream advection 
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Table 2 

Measured full widths at half max (FWHMs) for all regions. 





FWHM (arcsec 

) 




rriE ( 

:-) 


Region 

Band 1 

Band 2 

Band 3 

Band 4 

Band 5 

Bands 1-2 

Bands 2—3 

Bands 3—4 

Bands 4—5 


(0.7-1 keV) 

(1-1.7 keV) 

(2-3 keV) 

(3-4.5 keV) 

(4.5-7 keV) 

(1 keV) 

(2 keV) 

(3 keV) 

(4.5 keV) 

1 


8.80+2:1® 


7-40tS:l3 

' -0.42 


-0.47 ± 

0.06 

0.38 ± 0.13 

-0.70 ± 0.22 

2 


4-22tS:l| 

2 Qc+O-12 
^•'JD_0,09 

3-ootS:l2 

4 11+0.34 
^•-^-^-0.30 


-0.84 ± 0.08 

0.59 ± 0.16 

0.77 ± 0.23 

3 


2 47+0'0® 
-0.07 

1 -70-^0.09 
-^•"°-0.07 

^•-^^-0.11 

1 q9 + 0.10 

l-’J-^0,09 


-0.47 ± 

0.08 

0.41 ± 0.17 

-1.15 ± 0.22 

4 

0.00_Q 33 

2±.OO_0 08 

3-26ti;:“ 


3 20'*'^'^^ 

-0.83 ± 0.18 

-0.41 ± 

0.05 

0.31 ± 0.11 

-0.35 ± 0.17 

5 


^•^^-0.12 

3.06;“;(( 

Q 9^+0.15 

3 04 


-0.56 ± 

0.06 

0.15 ± 0.14 

-0.17 ± 0.19 

6 


9 QO"I“0.05 

OQ 

z.i?o_Q 09 

Z.UO_o 09 

9 91 -1-0.15 

^•^-'--0.14 

-0.19 ± 0.21 

0.36 ± 

0.06 

-0.92 ± 0.13 

0.18 ± 0.19 

7 

^•oy_o,i7 

ry Qn-|-0.05 

^•90_o,o5 

901 +0.08 

^•'Jl_0,08 

1 01 -l-O - 09 
-'-•°-'--0.07 

1 88+0 -11 
l.OO_0 08 

-0.39 ± 0.20 

-0.01 ± 

0.06 

-0.60 ± 0.14 

0.02 ± 0.17 

8 

r\ qQ+ 0'21 

Z.OO_o,20 

ry -1-0.08 

' ^-0,08 

ry QQ+0.10 

^•oo_0.09 


9 q'7 + 0.20 

-0.17 

0.43 ± 0.26 

-0.19 ± 

0.07 

-0.30 ± 0.15 

0.29 ± 0.22 

9 

2 16 + 0'24 

9 oi:; + 0'07 

Z.OO_o 06 

2 47+0'ii 
-0.11 

1 91+0-09 
-^•^-^-0.09 


0.24 ± 0.31 

0.07 ± 

0.07 

-0.63 ± 0.16 

0.34 ± 0.22 

10 

9 qq + 0.24 
Z.OO_o 23 

x.y3_Q 06 

1 70+0.09 

J-- ' 'J 0.08 

T.u»_0 08 

1 co+0.13 
-‘-•^°-0.12 

-0.50 ± 0.29 

-0.18 ± 

o.os 

-0.24 ± 0.18 

-0.02 ± 0.23 

11 


Q Oq-I-O. 15 

o--^9_0,13 

9 co+O.lO 
^•0^-0.13 

I.yu_0,i3 

'^•'^^-0.38 


-0.36 ± 

0.10 

-0.70 ± 0.22 

1.21 ± 0.37 

12 


3-86t“;l^ 

^•0-‘-_o.i3 

q r\9+0-22 

9 9q + 0.21 

'^•^'^-0.17 


-0.56 ± 

0.10 

0.36 ± 0.22 

-0.74 ± 0.27 

13 

2.85tm 

2-43ti::S® 

9 q/:?-|-0.08 
^•'30 o 05 

-L-y9_o.io 

-L-0^_0.14 

-0.45 ± 0.20 

-0.04 ± 

0.05 

-0.47 ± 0.13 

-0.15 ± 0.20 

14 

2-86tS;l^ 

^•^^-0.04 

9 9 Q+0.08 

^•^■3_0,07 

9 qO+O.lO 

Z.OO_o 08 

2 19+0-12 
'^•-^^-0.10 

-0.47 ± 0.17 

-0.12 ± 

0.06 

0.17 ± 0.12 

-0.20 ± 0.15 

15 

0 '71 +0.17 
' ^-0,16 

1 99+0'0® 

x.C7c»_q 04 

1 Qfi+O.OO 
l-oU_0 05 

00 

-4 

1 + 
0 0 
b b 

w -j 

1 c:9+0-09 

l-O^_0,08 

-0.85 ± 0.18 

-0.15 ± 

0.05 

0.09 ± 0.11 

-0.51 ± 0.16 

16 

1 q'7-1-0.14 

1-0'_0.13 

1 -70-1-0.04 
-*-• ‘'^-0.03 

1 r9 + 0.06 

-^•O^_0,05 

1 9r:+0-06 

-'-•^O_o,o4 

1 9q-l-0.08 

l-^9_o,06 

-0.22 ± 0.21 

-0.18 ± 

0.06 

-0.49 ± 0.13 

-0.04 ± 0.17 

17 

i.oo_Q 2^2 

-^•==^■^-0.05 

1 54+0'0® 
-^•O^_0,07 

1 45+°-°^ 
-L-^O_0.06 

2 05"*”*^'^® 

0.43 ± 0.22 

-0.31 ± 

0.07 

-0.16 ± 0.15 

0.86 ± 0.21 

18 



Q 1 Q+0'17 

d- J-O_o.i6 

z.yo_o 19 

-‘-•°^-0.16 


-0.49 ± 

0.09 

-0.17 ± 0.21 

-1.45 ± 0.32 

19 


Q Qfl-1-0.08 

^•oU-0.06 

9 98+0• 11 
^•^°-0.08 

2 10+0-12 

1.DU_o,i4 


-0.02 ± 

O.OS 

-0.13 ± 0.17 

-0.74 ± 0.27 

20 

4 81+"'^^ 

^•C>l_Q 31 

i-o^_0.03 

1 0-7+0.08 
-0.06 

1 =ifi+o-o^ 
-'-•O0_o.06 

2 14+0-23 
'^•-‘-^-0.23 

-2.68 ± 0.19 

0.02 ± 

0.07 

-0.44 ± 0.14 

0.77 ± 0.28 

Mean 

2.89 ± 0.35 

3.11 ± 0.37 

2.53 ± 0.23 

2.47 ± 0.30 

2.35 ± 0.23 

-0.46 ± 0.24 

-0.25 ± 

0.06 

-0.14 ± 0.10 

-0.09 ± 0.15 

Note. 

— Mean values computed for all regions; 

mean values are averages 

for region ttie 

values (i.< 

3., not 

computed from 

mean FWHMs) 


Errors on mean values are standard errors of the mean. Horizontal rules group individual regions into filaments. 


Table 3 

Best width-energy fit parameters, 

M = 1(2 = 1 

Table 4 

Best model fits for all regions, r ?2 derived from srcut fits, n 

Region 

Loss-limited 

Damped 


Regioi 


Loss-limited 

Damped 


Bo (a«G) 

x" Bo (mG) 


<2{j 

1 V2 Bq (fiG) 

Bo (mG) 

X^ 

<2{j 

1 

182 

35.1 

28 

24.3 

0.008 

1 

12.0 

256 

29.1 

19 

25.5 

0.008 

2 

312 

80.9 

25 

74.4 

0.003 

2 

11.8 

443 

99.4 

18 

72.0 

0.003 

3 

426 

21.5 

27 

23.9 

0.002 

3 

10.6 

599 

31.8 

19 

29.3 

0.002 

4 

284 

50.9 

28 

29.9 

0.004 

4 

12.0 

402 

36.2 

19 

34.5 

0.004 

5 

288 

19.3 

25 

14.5 

0.003 

5 

12.9 

419 

32.5 

19 

19.1 

0.004 

6 

410 

173.7 

68 

65.2 

0.004 

6 

11.9 

568 

89.7 

147 

61.9 

0.006 

7 

418 

50.9 

291 

11.6 

0.006 

7 

15.3 

617 

11.7 

473 

11.1 

0.009 

8 

388 

47.7 

138 

10.0 

0.005 

8 

8.4 

511 

15.1 

333 

10.0 

0.008 

9 

414 

65.1 

78 

14.2 

0.004 

9 

11.2 

573 

23.3 

66 

13.7 

0.005 

10 

466 

17.1 

29 

1.5 

0.002 

10 

9.4 

631 

1.4 

478 

2.0 

0.007 

11 

355 

12.0 

317 

15.6 

0.010 

11 

6.1 

452 

18.1 

21 

17.5 

0.003 

12 

317 

9.6 

26 

10.4 

0.003 

12 

3.7 

375 

12.6 

21 

9.7 

0.003 

13 

400 

52.7 

258 

10.4 

0.006 

13 

8.4 

529 

13.8 

368 

9.9 

0.008 

14 

383 

102.9 

60 

12.7 

0.004 

14 

11.9 

547 

28.6 

120 

12.4 

0.006 

15 

431 

70.3 

48 

20.0 

0.003 

15 

9.5 

598 

23.1 

426 

20.0 

0.007 

16 

493 

15.0 

434 

4.0 

0.006 

16 

8.4 

670 

6.3 

22 

5.8 

0.002 

17 

467 

61.6 

232 

29.0 

0.004 

17 

9.0 

637 

32.6 

32 

29.4 

0.003 

18 

283 

28.1 

24 

27.5 

0.003 

18 

14.6 

432 

72.4 

17 

16.6 

0.003 

19 

401 

36.9 

78 

8.3 

0.004 

19 

8.9 

542 

12.8 

209 

7.7 

0.006 

20 

463 

121.8 

58 

93.8 

0.003 

20 

9.6 

631 

96.0 

29 

94.2 

0.003 


Note. — Fits for Regions 1-3, 5, 11, 12, 18, and 19 have 3 degrees Note. — 772 values are computed from equation and held fixed 

of freedom; all others have 4. The choice of a best at value may be in model fits. All comments for Table apply to srcut fits as well, 

construed as removing one additional dof. Damped fits require at to 
be smaller than the FWHM at 2 keV in order to rule out effectively 
loss-limited fits with large a^. 
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Figure 5. Rim width fits as a function of energy for loss-limited and damped fits with /i = 1 and 772 = 1 fixed for all regions, from best 
fit parameters in Table Damped model predictions (e.g., Region 5) are not given at low energies if FWHMs cannot be calculated for 
model profiles (i.e., mooted intensity behind thin rim exceeds half-maximum of rim peak within model domain of ~20''). Ordinate (y) 
axis limits vary between subplots and are offset from the origin to better show model predictions and data variation. 


through a decreased effective velocity Vd — {dD/dx) in 
our transport equation: 


Vd - 


dD{x) 


dx 


= (14) 


requiring thinner rims as 772 and hence dD{x)/dx in¬ 
crease. Moreover, the effective velocity may impact rim 
widths even if advection dominated. In practice, if 772 


varies freely in damped fits, we observe that smaller val¬ 
ues of Of, permit and favor smaller best fit values of both 
Bo and 772 . 

Best fit Bq values are smallest for 772 approaching 0 
and /J. = 2 ; intuitively, 77 > 1 strengthens diffusion at en¬ 
ergies > 2 keV for fixed 772 , permitting a smaller best fit 
772 and hence smaller Bq- Fixing 772 = 1 (Tablehow¬ 
ever, ties the range of Bq values to the range of ^served 
rim widths. Our minimum loss-limited values of Bq 
are consistent with prior estimates of ~200-300 /rG for 
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Figure 6. Model predictions illustrate weak-field (Region 1) and 
strong-field damping (Region 16) for damped best fit parameters 
with p = 1 and r ;2 = 1 fixed. In X-ray energies (0.7—4.5 keV) model 
profiles are not strongly energy dependent, but weak-field damping 
profiles evolve and show no measurable FWHM at sufficiently low 
energies. Profiles are normalized to peak thin rim intensity, and 
shaded regions indicate damping lengthscale ah. Model parameters 
are given in Table 


advection-dominated transport (IVolk et al. 

2005 

Pari- 

zot et al.|2006 

Morlino & Gaprioli|2(J12p. It loss-limited, 


'I'ycho’s rims require strong magnetic Held amplification 
to ^lOOx typical galactic field values of ^2-3 /iG, versus 
the expected 4x amplification from a strong shock. 


Magnetic damping fits permit much smaller values of 
Bq, as expected. The minimum value of Bq is around 
20 /iG, which would require no field amplification beyond 
that from a strong shock. Fixing Bmin = 2 ^G instead 
of 5 /tG permits smaller fit Bq values, though fit values 
for both r ]2 and Bq still display considerable scatter. 

5. DISGUSSION 

5.1. Damping is poorly-constrained by X-ray 
width-energy dependence 

Width-energy dependence is not as sensitive a discrim¬ 
inant between damping and loss-limited models as may 
be intuitively expected. Both loss-limited and damped 
fits do not perfectly capture sharp drop-offs (especially 
between 0.7-1 keV and 1-2 keV), and the data show 
large scatter (Figure [^. In most regions, best fit width- 
energy curves have ttie —0.2. Only sharp decreases 
in rim width (ttie ~ —0.5) can disfavor damping as a 
control on rim widths (e.g.. Region 18). 

When width-energy dependence is weak, both loss- 
limited and damped models give best fits with similar 


profiles - amplified magnetic fields in both models cause 
rim intensities to drop sharply. When width-energy de¬ 
pendence is stronger, damping profiles yield energy de¬ 
pendent rims if magnetic fields are small (say, < 50 fxG) 
and the damping lengthscale is much smaller than rim 
FWHMs. Figure|^shows this contrasting behavior using 
best fit parameters for Regions 1 and 16. 

The best damped fit parameters for Regions 1-5, 10, 
12, and 18 predict a 1.375 GHz radio profile that does 
not drop below 50% of peak intensity within ^20" down¬ 
stream of the shock, which we treat as having no mea¬ 
surable FWHM. These regions are all best fit with field 
Bq < 40 pG. We refer to this model behavior as 
“weak-field” damping, associated with weak magnetic 
fields and stronger emission intensity immediately down¬ 
stream of the thin rim. As observation frequency de¬ 
creases, rim and trough contrast decreases, causing rim 
FWHMs to increase and eventually become unmeasur¬ 
able. Rim width-energy dependence strengthens dra¬ 
matically (toe —>■ —oo), permitting model fits to repli¬ 
cate strong width-energy dependence in observed rim 
widths. Weak-field fit parameters require energy losses 
downstream of the shock to produce energy dependence 
in a damped field, as discussed in Section 2.5.2| 

Damping model fits to all other regions predict consis¬ 
tently thin rims with measurable FWHMs at decreasing 
energy. The width-energy dependence parameterized by 
Toe trends towards zero at low energy, indicating that 
rim widths are comparatively energy independent. At 
higher energy, such rims narrow slightly with increas¬ 
ing energy (toe —0.2) to match the observed width- 
energy dependence in X-rays. The gradual increase in 
I Toe I with energy is expected in the damping model as 
advection and/or diffusion take control of rim widths at 
increasing energy, as discussed in Section [Z5| The best 
damped fits for these regions all have larger Bq values 
than in the “weak-field” damping case. We refer to these 
fits as giving rise to “strong-field” damping, associated 
with stronger magnetic fields, weaker emission intensity 
behind the thin rim, and clear rims with measurable 
FWHMs at low photon energies (at and below soft X- 
rays). 

We are unable to determine whether weak- or strong- 
field damping is descriptive of Tycho’s shock magnetic 
field, given the large values on our model fits. Sev¬ 
eral regions may be well-fit by “weak-field” and “strong- 
field” damping alike. But, the qualitative behavior of Bq 
is suggestive; very strong fields at the shock (> 100 pG) 
with magnetic field damping on a scale comparable to 
rim FWHMs appear incompatible with significant energy 
dependence. The distinction between weak- and strong- 
field damping may be equally well described as strong 
versus weak damping. The dichotomy reflects additional 
degeneracy between Of, and Bq - for a given rim width, 
increased damping length at requires increased Bq to re¬ 
produce the same width, and vice versa. 

Although model fits are poorly-constrained, we at¬ 
tempt to further explore how damping and synchrotron 
losses set rim widths for damped fits with finite Ob- 
Table compares rim widths at 2 keV to advection 
lengthscale Gd and damping lengthscale Ob from the loss- 
limited and damped fits of Table As discussed in 
Section [2.5[ loss-limited rim widths are qualitatively set 
by either damping or the dominant transport process as 
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Table 5 

Lengthscale analysis 


Region 

Measurements 

Loss-lim. fit 


Damped fit 


Vd 

(10® cm/s) 

w (2 keV) 
(%rs) 

^ad 

(%r.) 

lad/ldm 

C-) 

^ad 

(%r.) 

CLb 

(%rs) 

^ad / 

(-) 

w{2 keV)/a6 
(-) 

1 

1.30 

2.64 

0.61 

1.40 

10.26 

0.80 

12.83 

3.30 

2 

1.29 

0.99 

0.27 

1.39 

11.54 

0.30 

38.48 

3.28 

3 

1.29 

0.74 

0.17 

1.38 

10.80 

0.20 

54.00 

3.72 

4 

1.28 

1.36 

0.31 

1.38 

10.08 

0.40 

25.20 

3.40 

5 

1.28 

1.27 

0.30 

1.37 

12.07 

0.30 

40.24 

4.25 

6 

1.28 

1.24 

0.18 

1.37 

2.61 

0.40 

6.53 

3.10 

7 

1.27 

0.96 

0.17 

1.36 

0.29 

0.60 

0.49 

1.61 

8 

1.27 

0.99 

0.19 

1.36 

0.90 

0.50 

1.80 

1.98 

9 

1.26 

1.03 

0.17 

1.36 

2.09 

0.40 

5.23 

2.57 

10 

1.26 

0.73 

0.14 

1.35 

9.21 

0.20 

46.05 

3.66 

11 

1.25 

1.05 

0.21 

1.34 

0.25 

1.00 

0.25 

1.05 

12 

1.24 

1.09 

0.25 

1.33 

10.80 

0.30 

35.98 

3.62 

13 

1.23 

0.98 

0.18 

1.32 

0.34 

0.60 

0.57 

1.64 

14 

1.14 

0.93 

0.17 

1.23 

2.84 

0.40 

7.10 

2.32 

15 

1.14 

0.75 

0.15 

1.22 

3.95 

0.30 

13.16 

2.50 

16 

1.13 

0.63 

0.12 

1.21 

0.14 

0.60 

0.24 

1.06 

17 

1.12 

0.64 

0.13 

1.20 

0.36 

0.40 

0.91 

1.61 

18 

1.15 

1.32 

0.28 

1.23 

11.29 

0.30 

37.64 

4.42 

19 

1.18 

0.95 

0.17 

1.27 

1.99 

0.40 

4.97 

2.37 

20 

1.17 

0.78 

0.14 

1.26 

3.08 

0.30 

10.28 

2.59 


Note. — All lengthscales computed at fiducial energy 2 keV from best fits with r }2 — 1 and — 1. Ratio of /ad/^diff is same for loss-limited and 
damped fits and depends only on plasma velocity and observation energy, as ^ad/^diff is independent of Bq for ^ — 1. 


w ~ min (ob, max (/ad, ^diff))- At 2 keV with 772 = 1 and 
/i = 1 fixed, the ratio lad/ldiB is nearly constant; vari¬ 
ation (1.2-1.4) arises solely from azimuthal variation in 
shock velocity. If rims are loss-limited and diffusion is 
negligible, we anticipate rim widths w = 4.6/ad, where 
the factor 4.6 may be derived assuming spherica l sym¬ 
metry and an exponential synchrotron emissivity |Ballet 


20061. At 2 keV, we find that loss-limited rim widths 


w ^ 5/ad due to diffusion. Enforcing w{2 keV)/ab > 1 for 
damped fits still permits /ad/ofc < 1 as w{2 keV) < 5/ad, 
with damping decreasing w(2 keV) below the expected 
loss-limited width at a given Bq. We may impose a 
tighter or looser bound, but our results should not be 
greatly affected given large uncertainty in our fits and 
associated degeneracy between Bq and Ub- 
Our modeling can also reassess the significance of rim 
width-energy dependence in the remnant of SN 1006 pre¬ 
sented by |R14[ We fit averaged filament widths in SN 
1006 measured by |R14| to our damping model using the 
same procedure as for Tycho. We fix Rmin = 5 mG, al¬ 
though best fit magnetic fields for SN 1006 are smaller 
than those of Tycho due to the much wider filaments 
of SN 1006. A more extensive search of parameter space 
produces acceptable damping (hyb rid) models with more 
rapid shrinkage than found in |R14[ Damped fits are com¬ 
parable to or better than loss-limited fits in 3 of 5 fila¬ 
ments in SN 1006. Fits to two filaments with strong en¬ 
ergy dependence (ms ~ —0.5) favor a loss-limited model 
with sub-Bohm diffusion (772 ^ 1), though sub-Bohm 
diffusion may be an unphysical result pointing to over¬ 
simplifications in the model. The width-energy depen¬ 
dence in SN 1006 (ttte ~ —0.3 to —0.5) is overall slightly 
stronger than in Tycho, and the best damped fits for SN 
1006 all fall into the “weak-field” case. The best fit Bq 
values are less than 40 /iG in the damped model, com¬ 
pared to ~ 100-200 fiG in the loss-limited model. I f thin 
radio rims in SN 1006 (Reynolds fc Gilmore||1986) indi¬ 


cate magnetic field damping on lengthscales comparable 
to X-ray filament widths, some magnetic damping is not 
incompatible with rim width narrowing. A future study 
of SN 1006 with recent multi-frequency Karl G. Jansky 
Very Large Array data (PI: D. Green) may verify whether 
radio and X-ray rims can be jointly described by an am¬ 
plified and subsequently damped magnetic field. 


5.2. Model assumptions and correctness 

Our models adopted a distance d to Tycho of 3 kpc, 
but es timates for Tycho’ s distance range between 2.3- 
4 kpc (|Hayato et al.(2010|). As a larger remnant distance 
would increase both physical filament widths and shock 
velocity estimates from proper motion, we may derive 
expected scalings for modeled 772 and Bq as functions of 
assumed distance d. The advective lengthscale, equa¬ 


tion (101, may be rearranged to obtain: 


Bq = (3.17 TiG) 

^ad 


Vd 

10® cm/s 
- 2/3 


0.01 kpc 


2/3 


hv 

1 keV 


-1/3 


or, more simply, 

Bq CX {Vdf^^ (/ad)"^^^ 


(15) 


(16) 


Both /ad and Vd scale linearly with remnant distance d 
and thus their effects cancel in determining the mag¬ 
netic field. If diffusion is th e primary control on filament 
lengthscales, equation (fTTl) yields: 


772 oc (/diff)' 


(17) 


Model fits with varying distance obey both scalings, 
772 oc and Bq constant. When comparing model fits 
with remnant distances of 3 kpc and 4 kpc, the devia¬ 
tion from the idealized scaling is < 1% for Bq and ^1- 
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5% for r] 2 . As varying remnant distance d leaves width- 
energy scaling toe invariant, the relative contributions of 
/ad and /diff should also be invariant. Then both length- 
scales should scale simultaneously with d, yielding the 
observed behavior. In the damped model, a larger dis¬ 
tance d will require larger physical damping lengths from 
magnetic turbulence. But, fitted values of ab should re¬ 
main unchanged as we report at in units of shock radius 
rs- 

The exponential variation in D{x) (equation may 
not be physically reasonable; as noted above, rhe as¬ 
sumption D(x) oc 1/B^{x) gives rise to sharp gradients 
in diffusion coefficient and even inverts the effect of 772 on 
our model profiles (where larger 772 can cause rim widths 
to narrow). Nevertheless, the modeled behavior appears 
physically reasonable. Only X-ray profiles are impacted 
by the assumption on D{x) as radio profiles assume no 
diffusion. And, model behavior driven by advection and 
magnetic damping - namely, rim width-energy depen¬ 
dence in a damped magnetic field - occurs for D = 0, and 
cannot be an artifact of our assumption that D{x)B‘^{x) 
is constant. 

We emphasize that additional counts from averaging 
measurements or selecting larger regions will likely not 
improve our ability to constrain Bq and 772 from width- 
energy modeling. This is easily seen from Figurel^as well 
as large values in Tables |^that reflect relativdy tight 
errors in our FWHM profile measurements (Table [^. 


5.3. Other potential constraints 

Rim width-energy dependence is a morphological man¬ 
ifestation of spectral s oftening downstream of the for¬ 
ward shock; previously, Cassam-Chenai et al. (2007) also 
sought to distinguish loss-limited and damped rims with 
a careful spectral study and 1-D hydrodynamical model. 
Although knowledge of radial spectral variation, in prin¬ 
ciple, fully determines rim profiles, there is not a clear 
relationship between observed spectral variation and rim 
widths alone. Our model, for example, replicates ob¬ 
served rim width variation but underpredicts the ob- 
se rved spectral variation . The more sophisticated model 
of Cassam-Chenai et al. (2007) similarly has trouble re- 
producing the observed radial gradient in spectral pho¬ 
ton index. Photon indices predicted in our model, for 
integrated model rim spectra similar to those of Figure 
and Table are also somewhat ill-constrained. Enforc- 
in g limits on acceptable mode l photon indices (e.g., Fig. 
16 Cassam-Chenai et al. 2007) could help constrain shock 
parameters or identify discrepancies between model as¬ 
sumptions and measurements. 

Estimates of magnetic field strength that d epend upon 


global models for Tycho’s evolution (e.g., Morlino & 
Caprioli 2012) may not help constrain our rim model 


results. In such models, assumptions on particle accel¬ 
eration and magnetic field evolution downstream of the 
shock (e.g., further damping or amplification at the re¬ 
verse shock and contact discontinuity) would affect the 
spatially integrated spectrum of Tycho. We have focused 
on effects immediately behind the forward shock, making 
as few assumptions as possible about particle accelera¬ 
tion and magnetic fields throughout the remnant. It may 
not make sense to extrapolate our results to estimate 
magnetic fields throughout the remnant or, similarly, to 
use global models to constrain our results given all as¬ 



Figure 7. Radio image of Tycho’s SNR at 1.375 GHz with linear 
scaling. Extraction regions (green) for joint radio and X-ray pro¬ 
file analysis overlay region selections for X-ray rim width analysis 
(Figure m . 


sumptions involved. 


6. JOINT RADIO AND X-RAY MODELING 

Thin radio rims spatially coincident with X-ray syn¬ 
chrotron rims may help constrain magnetic field damp¬ 
ing. Radio synchrotron emission in the remnant interior 
may arise from not only recently shock-accelerated elec¬ 
trons but also long-lived electrons interacting with a tur¬ 
bulent field inside the remnant, as ^GeV electrons have 
cooling times of order 10^-10^ yr in 10-100 pG magnetic 
fields. But simple models for synchrotron emissivity as 
a function of density in Sedov and pre-Sedov dynami¬ 
cal stages predict radio synchrotron profiles rising grad- 
ually to broad maxima well inside the shoc k radius (e.g. 
Reynolds & Chevalierll981||Reynolds|1988|). Our steady- 
state planar transport model predicts monotonically in¬ 
creasing emission downstream of the forward shock in 
a loss-limited model, until sphericity becomes important 
and our model is inapplicable. On scales of a few percent 
of the remnant radius, only magnetic field damping can 
cause emission to decrease. 

The idea of mod eling X-ray and radio p r ofiles jointly 
was pioneered by |Cassam-Ghenai et al.| (|2007|), who 
found that sharp radio rims in Tycho were not repro¬ 
duced by a loss-limited model. Gassam-Chenai' et al.| 
(2007) considered loss-limited and damping profiles con¬ 


sistent with physical constraints applied in a hydrody¬ 
namical model. Here, we neglect physical constraints 
and estimate a damping length necessary to generate ra¬ 
dio rims with shape similar to observations. 

We extract radio profiles from a 1.375 GHz image of 
Tycho taken with the Very Large Array (VLA) in A con¬ 


figuration in March 1994 (PI: D. Moffett); see Reynoso 
et al. (1997) for a detailed presentation. The halt-power 
beam width of ~1.5" just resolves thin radio rims and 
structure near the forward shock; the image is sampled 
at 0.5". We also extracted 4-7 keV X-ray profiles in 
all regions from the previous archival Chandra observa¬ 
tion to jointly model radio and X-ray profiles, permit- 
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Table 6 

Fit parameters from profile shape comparison 


Region 

Bo 

(mG) 

Ub 

(-) 

Region 

Bo 

(mG) 

Qfb 

(-) 

A 

50 

0.020 

I 

250 

0.020 

B 

200 

0.050 

j 

300 

oo 

c 

15 

0.010 

K 

400 

0.010 

D 

100 

0.050 

L 

250 

oo 

E 

120 

0.030 

M 

200 

oo 

F 

300 

0.025 

N 

800 

0.010 

G 

250 

0.020 

o 

200 

0.005 

H 

300 

0.020 

p 

150 

0.012 


Note. — Damping lengths of oo indicate that a loss-limited fit is 
favored {ay > 10% of shock radius Vs). 


ting somewhat firmer discrimination of plausible model 
parameters. Figure shows the extraction regions over¬ 
laying the radio image and the previous X-ray profile 
regions of Figure 

We compute model radio and X-ray measured profiles 
from the transport model of equation Q for varying Bq 
and Qb, similar to those shown in Figure The pa¬ 
rameters = 5 /rG, 772 = 1, and p, = 1 are held 

fixed. Diffusion, in particular, is negligible for modeled 
radio emission as particle energies are 3 orders of mag¬ 
nitude lower than in X-ray. Neglecting diffusion also cir¬ 
cumvents the unphysical assumption that D{x)B'^{x) is 
constant, which was invoked to obtain Green’s function 
solutions to equation Q. 

We align each set of model profiles at some Bq and 
ab to measurements by eye, varying relative amplitudes 
and translations in radio and X-ray independently to 
best match the measured profiles. Although “fitting” 
by eye does not quantitatively bound model parame¬ 
ters, we can find plausible values of Bq and ab and es¬ 
timate the importance of magnetic damping throughout 
the remnant. More involved nonlinear fitting may be un¬ 
reliable as we cannot constrain spatially heterogeneous 
radio emission within the remnant. Moreover, our trans¬ 
port model neglects self-similar downstream evolution of 
shocked plasma (decaying velocity, density, and magnetic 
field) and is inaccurate further downstream than ^10% 
of the shock radius. 

The joint radio and X-ray profile modeling contrasts 
strongly with our previous width-energy fitting from X- 
ray measurements alone, which neglected profile shape in 
favor of more robust FWHM measurements. Manually 
fitting profiles allows us to consider radio and X-ray fila¬ 
ments that do not have well-defined FWHMs, especially 
as radio rims do not fall below 50% of the peak emission. 
We can also use profiles from regions not previously con¬ 
sidered due to the lack of an X-ray FWHM, especially in 
softer (0.7-4 keV) X-rays. 

We identify three classes of radio profiles: thin rims 
with downstream troughs, plateaus, and continuous rises. 
Regions B, G, and D (southern limb) show plateaus in 
radio emission. Regions J, L, and M (around NW) show 
continuous rises in emission. All other regions have a 
radio rim within 15" of the forward shock, where the for¬ 
ward shock in radio is assumed at zero intensity. Figure!^ 
shows extracted radio and nonthermal X-ray (4-7 keV) 
profiles for three well-fit regions with our best manually 
selected model profiles, illustrating each of the three ra¬ 
dio profile types observed. 
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Figure 8. Measured radio and X-ray profiles plotted with model 
profiles for varying a^, and Bq in each region, showing typical pa¬ 
rameters (and ranges) required to reproduce radio and X-ray rims 
simultaneously in our model. Solid black curves in all regions plot 
our manually chosen best model profiles. Profiles are chosen to 
show varying radio rim morphology, including plateaus (B), thin 
rims with troughs (K), and continuous rises (M); these profiles 
show some of the best agreement of our selected regions, but cf. 
Figure Profile radial coordinates are shifted arbitrarily to aid 
visual comparison. Negative radio intensity is unphysical and as¬ 
sociated with deconvolution of raw VLA visibilities. 


Our model requires damping length ab < 0.1 to pro¬ 
duce a plateau or thin rim in radio emission. For regions 
with thin radio rims, the best manually selected profiles 
have Bq between 50-400 /rG, neglecting only Region N 
which could not be modeled simultaneously in both X- 
ray and radio (Figure^. The damping length ab ranges 
between 0.01-0.03, or z-7". We list estimated best fit 
parameters in Table 

The radio plateaus in regions B, G, D are compatible 
with damping lengths between 0.01-0.05. Gontinuous 
rises in radio emission are best modeled with ab ^ 0.1 and 
i?o^200-300 TiG. Although damping lengths ab > 0.1 
are not physically meaningful well beyond the shock, 
these large values of ab yield practically constant mag¬ 
netic field near the shock. Our results are also insen¬ 
sitive to the assumed value of Bmin; only model pro¬ 
files with small magnetic fields {Bq < 100 fxG) and/or 
small damping length {ab < 0.01) are affected if we take 
Bmin = 2 fJ,G rather than 5 /rG. 

Tycho’s shock structure is more complex than assumed 
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Figure 9. Extracted profiles poorly reproduced by our model, 
compared to Figure Region I shows irregularly shaped radio 
rim. Region N contains two superposed filaments that cannot be 
modeled by a single rim in both X-ray and radio; the narrow X- 
ray rim requires atypically strong magnetic fields (~ 800 ^G), and 
the emission plateau behind the radio rim requires small damping 
lengths (ctf, 0.005). 


by our transport model. Emission towards the remnant 
interior clearly shows spatial structure (Figurej^. Figure 
shows two regions that were poorly described by our 
model. A majority of our regions have irregular rims; 
in at least 2-3 regions, this may be attributed to pro¬ 
jection of multiple filaments. Others (e.g., Region K) 
show rims with slopes that cannot be matched by our 
models, whether too steep or shallow. Shape mismatch 
may be attributed in part to point-spread mismatch, dif¬ 
fusion (e.g., 772 ^ 1), shock precursors, or other effects. 
Nevertheless, the conclusion that thin radio rims require 
magnetic damping is supported by more sophisticated 
modeling. As we have noted, hydrodynamic models with 
diffusive shock acceleration (Cassam-Chenai et al.|[20R7 


Slane et al.|2014 ) also cannot produce radio prohles with 


narrow rims m a purely advected magnetic field. It is 
also not clear that radio structure is due to a global ra¬ 
dial va riation in magnetic field strength. E.g., |Slane et al.| 
(2014) suggest that Rayleigh-Taylor fingers between the 
forward shock and contact discontinuity could locally 
confine radio-emitting electrons, creating observed radio 
rim structure. But, if radio rims are due strictly to mag¬ 
netic field effects, our conclusions on magnetic field drop¬ 
off are independent of the causative physical mechanism 
so long as they yield magnetic field fall-off over 1% of 
remnant radius. 

We may also use spatially overlapping radio and X- 
ray region selections (Figure to attempt to constrain 
X-rw width-energy fits. All Xu-ay region selections (Fig¬ 
ure [R are associated with radio rims, except for Regions 
11 and 12 underlying Region J in radio. X-ray emission 
at regions B, C, D (radio plateaus) and L, M (radio rises) 
showed rims, but their widths either could not be mea¬ 
sured, or could only be measured in the 4-7 keV band. 


If radio rims require damping at ^1% remnant radius, 
width-energy fits with “weak-field” damping lengthscales 
of ~0.5% may be disfavored. Conversely, the radio and 
X-ray filament of Regions 11,12 and J is the best remain¬ 
ing loss-limited rim candidate, although width-energy fits 
are equivocal towards damped and loss-limited models. 


7. CONCLUSIONS 

We measured the widths of several thin synchrotron 
filaments around Tycho’s supernova remnant and found 
moderate narrowing of rim widths throughout the rem- 
nant, corrob orating rim narrowing observed by |Ressler| 
et al. (2014) in the remnant of SN 1006. We confirmed 
that selected filaments are dominated by nonthermal 
emission and have clearly measurable full widths at half 
maximum in 4-5 energy bands. Both X-ray width-energy 
fits and joint radio/X-ray profile modeling require mag¬ 
netic fields > 20 /rG even with magnetic damping. 

A steady-state particle transport model with constant 
magnetic field gives diffusion coefficients and magnetic 
field strengths broadly consistent with prior estimates 
from rim widths (e.g., Parizot et al. 2006 Rettig fc 


Pohl|201^ and r adio and gamma ray m easurenients (Ac- 
ciari et ah] 2011 [Morlino fc Caprioli 12012 ). The same 
model with a damped magnetic held is equally capa¬ 
ble of describing our measured data. At weak energy 
dependence the two models are indistinguishable and 
magnetic damping fits favor moderately amplified mag¬ 
netic fields beyond simple compression, but lower than 
for loss-limited models. At moderate energy dependence 
(toe —0.3), the damping model permits weak magnetic 
fields and short damping lengths (< 1% remnant radius) 
to reproduce energy dependence, but we still cannot fa¬ 
vor either damped or loss-limited rims due to fit uncer¬ 
tainty. The distinction between loss-limited and damped 
models is somewhat artificial; only for large differences 
between the various transport lengthscales, a frequency- 
dependent occurrence, is one or the other mechanism 
clearly dominant. 

Thin radio synchrotron rims, however, are not re¬ 


produ ced in loss-limited models (Cassam-Chenai et al. 
2007). Assuming shocked electrons account for most ra¬ 


dio emission immediately downstream of Tycho’s forward 
shock, we jointly model radio and X-ray profiles and find 
that damping lengths of 1-5% of the shock radius are 
required throughout most of the remnant; only a few 
(3/16) selected regions are plausibly consistent with a 
constant advected magnetic field. Typical magnetic field 
strengths range between 50-400 /iC. Although we can¬ 
not bound damping lengths and fields from our qualita¬ 
tive profile comparisons, our results are physically rea¬ 
sonable and are likely good to order-of-magnitude. If 
damping lengths inferred from radio rims are correct, 
“weak-field” damping is disfavored in explaining X-ray 
rim width-energy dependence, and damped rims require 
magnetic field amplification to ~100 /rC or more in Ty¬ 
cho. 
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The scientific results reported in this article are based 
on data obtained from the Chandra Data Archive. This 
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research made extensive use of NASA’s Astrophysics 
Data System. This research also made use of APLpy, 
an open-source plotting package for Python hosted at 
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